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Abstract 

In the present work we study the effects of an unparticle IA as the possible source of 
missing energy in the decay B — > K{K*) + missing energy. We find that the dependence of 
the differential branching ratio on the K(K*)-meson , s energy in the presence of the vector 
unparticle operators is very distinctive from that of the SM. Moreover, in using the existing 
upper bound on B — > K(K*) + missing energy decays, we have been able to put more 
stringent constraints on the parameters of unparticle stuff. 
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1 Introduction 



Flavour Changing Neutral Current (FCNC) processes are not only powerful tests of the Standard 
Model (SM) but also provide very stringent tests for any physics beyond it. The smallness of 
FCNC processes in the SM is attributed to the fact that these processes are generated at loop 
level and are further suppressed by the CKM factors. Due to their smallness within the SM 
these processes can also be very sensitive to any new physics beyond the SM. Amongst the many 
FCNC decays involving B and if-mesons the decays of the form b — > s + missing energy have 
been the focus of much investigation at the B factories Belle and Babar. 

Of particular interest, in the SM, is the decay b — > suu, as it has the theoretical advantage of 
uncertainties much smaller than those of other decays, due to the absence of a photonic penguin 
contribution and hadronic long distance effects. However, in spite these theoretical advantages, it 
might be very difficult to measure the inclusive mode B — > X s uu, as it requires a construction of 
all the X,'s. Therefore the rare B — > K(K*)uu decays play a special role, both from experimental 
and theoretical points of view. Also the branching fractions of the 5-meson decays are quite 
large, with theoretical estimates of Br(B — > K*vv) ~ 1CT 5 and Br(B — > Kvv) ~ 1CT 6 [1]. 
These processes, based on b — > suu, are very sensitive to non-standard Z models and have been 
extensively studied in the literature [2-4]. 

As such, any new physics model which can provide a relatively light new source of missing 
energy can potentially enhance the observed rates of B — > K(K*) +missvs\g energy (B — > K(K*)+ 
ft), where many models have been proposed which provide such low mass candidates (which can 
contribute to b — > smissing energy). Note that in reference [3] the phenomenology of such 
low mass scalars was explored. Such studies have also been done in the context of large extra 
dimension models [5] and leptophobic Z' models [1,2]. One such model, which has excited much 
interest recently, is that of Unparticles, as proposed by H. Georgi [6]. In this model we assume 
that at a very high energy our theory contains both the fields of the SM and the fields of a theory 
with a nontrivial IR fixed point, which he called the Banks-Zaks (BZ) fields [7]. In his model 
these two sets interacted through the exchange of particles with a large mass scale M u , where 
below this scale there were nonrenormablizable couplings between the SM fields and the BZ 
fields suppressed by powers of M u . The renormalizable couplings of the BZ fields then produced 
dimensional transmutation, and the scale-invariant unparticle fields emerged below an energy 
scale A u . 

In the effective theory below k u the BZ operators matched onto the unparticle operators, 
and the nonrenormaliable interactions matched onto a new set of interactions between the SM 
and unparticle fields. The end result was a collection of unparticle stuff with scale dimension 
dui which looked like a non-integral number du of invisible massless particles, whose production 
might be detectable in missing energy and momentum distributions [6]. 

Recently there has been a lot of interest in unparticle physics [6, 8-16], where the signatures of 
unparticles have been discussed at colliders [8, 10, 15], in Lepton Flavor Violating (LFV) processes 
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[13], cosmology and astrophysics [16], and low energy processes [9, 11, 12]. 

In the present work we study the B — > K(K*)+J£ decay in unparticle theory, where this work 
is organized as follows: In section 2 we calculate the various contributions from both the SM and 
unparticle theory to the above-mentioned decays. Section 3 contains our numerical analysis and 
conclusions. 



2 Differential Decay Widths 

In the SM the decay mode B — > K(K*) + ft is described by the decay B — > K(K*)uu. As 
was noted earlier, unparticles can also contribute to these decays. Hence a comparison of the 
signatures of the two decay modes B — > K(K*)uu and B — > K(K*)U is required. 

In the SM the decay B — > K(K*)vi> is described by the quark level process b — > si/P through 
the effective Hamiltonian: 

W = ^Wio *7#. t 1 " 7s ) (1 - 7s) ^ , (i) 

where 

C = 4^- , (2) 
and the X(a; t ) is the usual Inami-Lim function, given as: 

with x t = m t /m 2 v . 

Similarly, the unparticle transition at quark level can be described by b — > sZY, where we shall 
consider the following operators: 

Scalar unparticle operators =>■ C s -^- s^/^b d^Oy + Cp-^- ^7^75^ <9 M O w , 

Vector unparticle operators =>- Cy d _ x s^^b Oy + Ca d _ x sj^^b Oy . (4) 

Before proceeding with our analysis note that we shall write the propagator for the scalar 
unparticle field as [8, 10]: 

/ d 4 xe iP *{0\TO u (x)O u (0)\0) = ^ J^^ {-P 2 ) du ~ 2 , (5) 

where 

16tt 5 / 2 T{d u + 1/2) 
du (2ir)^T(d u -l)T(2d u ) ' 
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2.1 The Standard Model 

Using the SM effective Hamiltonian for the quark level process b — > suu, as given in equation (1), 
we can calculate the differential decay width of B — > K{K*)vv (using the form factor definitions 
for the B — > K transition as given in appendix A.l). 

After taking into account the three species of SM neutrinos, we evaluate the differential decay 
width as a function of if-meson energy (E K ) as: 



dT SM _ G\a 2 
dE K ~ 2 7 7r 5 m| 



V ts V t l\ 2 \C 10 \ 2 fl(q 2 ) \ 3/2 (m 2 B ,m 2 K ,q 2 ) , 



(6) 



where X(m B , m 2 K , q 2 ) = m B + mj < + q 4: — 2m 2 B q 2 — 2m 2 K q 2 — 2m 2 K m 2 B , and q 2 = m 2 B + m 2 K — 2m B EK- 
Similarly, for the B — > K* case, using the definition of form factors for B — > K* transitions 
as given in appendix A. 2, the differential decay rate in the SM can be calculated as: 



dY SM _ G 2 F a 2 
dE K , ~ 2% 5 m| 



V ts V: h \ 2 \V 2 \C w \ 2 U\q 



V 2 



+ 



(m B + m K *) 2 m 2 K 



A 2 



A 2 



(m B + m K *) 2 



+ (m B + m K *Y(X + \2m 2 K ,q 2 )A{ - 2X(m 2 B - m\* - q 2 )Re(A{A 2 ) 



where A = m% + m K * + q — 2m B q — 2m K *q — 2m Kt m B , and q = m B + m K * — 2m B Ex 



(7) 



2.2 The Scalar Unparticle Operator 

As listed earlier, the following scalar operators can contribute to the B — > K(K*) U decay: 



HA 



HA 



where we have defined our form factors in appendix A. As such, the matrix element for the 
process B(p) — > K{jpf) + U(q) can be written as: 



Ate 



C, 



f+(m 2 B -m 2 K ) + f-q 2 



(9) 



The decay rate for B(p) — > K(p')U(q) can now be evaluated to be: 
dT su 1 A d , 



dE K 8ir 2 m B A 2 f u 



& \Cs\ 2 \lE 2 K -m 2 K (m% + m 2 K - 2m B E K ) 
u 



d u -2 



X 



f+(m B - m 2 K ) + f-(m 2 B + 2m 2 K - 2m B E K ) 



1 2 



(10) 



For the B — > K* transition our calculation proceeds along the same lines, where the matrix 
element for B(p) — > K*(p')U(q) can be written as: 



iC 



M s = -^J(e-g) { ( m B + m K *)A 1 - (m B - m K *)A 2 - 2m K * (A 3 - A )} O u 



(11) 
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and the differential decay rate as: 

dT SU m B A du 2 2 / 2 2 \ 3 / 2 / 2,2 r> 7-i \ d «- 2 / 10 n 

Zz^T = 2^ ' P ' ^* ~ mK * J \ mB + mR * ~ 2m B E K* ) • (12) 

As can seen from the above expressions the scalar unparticle contribution to the decay rate 
for B — > KU and B — > K*U will depend upon C s and Cp respectively. This shall allow us to 
place constraints upon C s and C P from these two different decay modes. This issue shall be 
re-visited in the final section of this paper. 



2.3 The Vector Unparticle Operator 

Along similar lines as followed in the previous subsection, we shall now make use of the vector 
unparticle operators: 

Cyp^r s^b C>£ + Caj^—t 57^756 0£ = s7m (<V + <?a7s) & 0£ , 

14 14 

and the form factors of appendix A, to calculate the matrix element for B(p) — > K(p')U(q): 

' ' r /+(p + p')M + /-(p-p')J . (13) 



.M V = — j- — f Cy 



And as such we calculate the differential decay rate as 

dT vu 1 A dll ,„,„„.,/,. , . ,. v /«-2 



dEft OTY~'fllB 



K~™>K 

x{-(m 2 B + m 2 K + 2m B E K ) + 2 T ~ P J ■ ( 14 ) 

I (mfj + mjf - 2m B E K ) J 

For the B K* case the matrix element for — > K*{p')U(q) is: 

+ "^1 ( ^ W^W 7 )) 0£ . (15) 

A^ w \m B + my. p r /J 

And therefore the differential decay rate will be: 

jpvw 1 , 4 r t/2 



8ttW ' (A^- 1 ) 1 V K V( mB + mi ,,) 

1 



2 



m 2 K ,(m B + m K *) 2 q 



[m B + mx*) 4 (377ix* + 2m B m 2 K , — 6m B m 2 K «EK* + m 2 B E 2 K *)A\ 



+4m B (E 2 K * - m 2 K ,) 2 A 2 2 + A(m B + ffl]f.)Vi!^« - m 2 K *)(m 2 K * - E 2 K Am 2 B A x A 2 



(16) 



To obtain the total decay width for B — > we must integrate over E K in the range 
m K < E K < (fflg + m 2 K )/2m B , whereas to obtain the total decay width for B — > we must 
integrate over E^* in the range m K , < E K * < {m 2 B + m 2 K A /2m B . 
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Figure 1: The differential branching ratio for: (a) Left panel: B — > K + E as a function of the 
hadronic energy (Ek). (b) Right panel: B — > K* + E as a function of the hadronic energy (Ex* )■ 
The other parameters are d u = 1.9, A u = lOOOGeF, C P = C s = 2 x 1(T 3 and C y = C A = 10 -5 . 



3 Numerical Results and Conclusions 

The total contribution to B — > K(K*) + |£ can be written as: 

T = T SM + T u , (17) 

where the T SM and T w are the SM and unparticle contributions to the B — > K{K*) + ^ decay. 
And we should note that in the SM the missing energy in the final state is attributed to the 
presence of neutrinos. Hence the SM contribution to this process is given by B — > K(K*)vi>. In 
the present case this signature can be mimicked by the process B — > K(K*) U, where we shall 
now try to estimate the bounds on the unparticles from the experimental constraints on missing 
energy signatures, as given by the 5-factories BELLE and BaBar [17, 18]: 

Br(B -> Kvv) < 1.4 x 1(T 5 , 
Br(B -> K*vv) < 1.4 x 10~ 4 . 

It is important to note that the SM process B — ► K{K*)vv provides a unique energy distri- 
bution spectrum of final state hadrons (K/K* in our case). Presently the experimental limits on 
the branching ratio of these processes are about one order below the respective SM expectation 
values. However, these processes are expected to be measured at future SuperB factories. As 
such, we presently only have an upper limit on the branching ratio of these processes, where to 
estimate the constraints on the unparticle properties. 

Note that H. Georgi, in his first paper on unparticles, tried to emphasize that unparticles 
behave as a non- integral number of particles [6]. He further went on to analyze the distribution 
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Figure 2: The branching ratio for B — * K + ft as a function of dy for various values of An. The 
left panel is for the contribution from the scalar operator, and the right panel is for the vector 
operator. The other parameters are Cs = 2 x and Cy = 10~ 5 . 



of the w-quark in the decay t — ► ulA. It was argued that the peculiar shape of the distributions 
of E u (the energy of the u-quark) may allow us to discover unparticles experimentally. As such, 
we have attempted to extend this same analogy to the process presently under consideration. 

Finally, before presenting our numerical results, note that the future SuperB factories will be 
measuring the process B — > K(K*) + ft by analyzing the spectra of the final state hadron. In 
doing this measurement at B-factories a cut for high momentum on the hadron is imposed, in 
order to suppress the background. Recall that unparticles would give us an unique distribution 
for the high energy hadron in the final state, such that in future 5-factories one will be able to 
distinguish the presence of a scale invariant sector (or unparticles) by observing the spectrum of 
final state hadrons in B — >• K(K*) + ft. 

With this idea in mind we have tried to plot the differential decay width of B — > K(K*) + 
ft as a function of Ek{Ek*) in figure (1). As we can see from these figures the unparticle 
operators (especially the vector operators) give us a very distinctive distribution for the final 
state hadron's energy. The distribution of the unparticle contribution is strikingly different when 
we include a vector operator for a highly energetic final state hadron. As such, unparticle stuff 
can give a distinctly different signature from the SM in this regime, which it should be noted is 
experimentally more favorable at future SuperB factories. 

In the next set of figures, figure (2), we have tried to analyze the constraints on the unparticle's 
scaling dimensions (dy) from different values of the cut-off scale Ay. In these plots we have used 
some specific values of the effective couplings Cs, Cp, Cy and Ca- As we can see from these figures 
the branching ratio is very sensitive to the scale dimension du and Au- In figure (3) we have 
shown the same plots for B — > K* + ft. From these two figures we can observe that the vector 
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operators are more strongly constrained as compared to scalar operators. The second feature is 
that B — > K + ft provides better constraints than the B — > K* + ft decay. 
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Figure 3: The branching ratio for B — > K* + ft as a function of du for various values of A^. The 
left panel is for the contribution from the scalar operator, and the right panel is for the vector 
operator. The other parameters are Cp = 2 x 10~ 3 and Cy = Ca = 10~ 5 . 



We have next tried to estimate the limits on the allowed values of the effective couplings, Cs, 
Cp, Cy and Ca, from the present experimental limits on the branching ratio of B — > K(K*) + ft. 
Therefore, in figure (4) we have shown the dependence of the branching ratio of B — > K + ft on 
Cs and Cy- As we can see from the expressions of the differential decay rate for B — > K + ft, 
given in the previous section, if we consider the scalar (vector) operators, then the rate for this 
process is only dependent on Cs (Cy). 

Finally, in figure (5) we have shown the dependence of the branching ratio of B — > K* + ft 
on the effective vertices. If we consider scalar operators then the rate of this process is only 
dependent upon Cp, whereas if we consider the vector operators then the rate can depend upon 
both Cy and Ca- 

To re-emphasize these last few points: 

• B — > K+ scalar unparticle operator shall constrain the parameter Cs, 

• _£?—> K*+ scalar unparticle operator shall constrain Cp, 

• B — > K+ vector unparticle operator will constrain only Cy, 

• whilst B — > K*+ vector unparticle operator will constrain both Cy and Ca- 

To conclude, in this work we have analyzed the effects of unparticles on the missing energy 
signatures of rare 5-decays. We have tried to argue that B — > K(K*) + ft provides very useful 
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Figure 4: The branching ratio for B — > K + ft as a function of Cs (left panel) and Cy (right 
panel). The cutoff scale has been taken to be k u = 1000 GeV. 



constraints on the parameters of the model, where we have considered four operators, namely 
the scalar, pseudo-scalar, vector and axial vector operators. Both the modes B — > K + ft and 
B — > K* + ft are different functions of these four operators, and hence provide independent 
constraints on the parameter space of the model. The set of operators we have considered, in 
principal, also contributes to meson anti-meson mixing, namely, K — K, Bd — Bd, B s — B s and 
D — D mixing. Detailed analyses of these, within the context of unparticle physics, has been 
done in reference [12]. Finally, note that the constraints provided by B — » K(K*) + ft in some 
regions can be much stronger than the ones provided by meson anti-meson mixing. 

Acknowledgement 

The work of NG was supported by JSPS under grant no. P-06043. NG would also like to thank 
Yasuhiro Okada and Sukanta Dutta for the discussions he had with them. We would also like 
to thank Steven Robertson and Kai-Feng Chen for their comments regarding missing energy 
signatures at S-factories. 

A The Form Factors 

A.l The form factors for the B — > K transition 

The form factors for the B — > K transition can be written as [19]: 

(K^)\s llM b\B{p)) = (p+pV + + ? m /_ , 
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Figure 5: The branching ratio for B — > K* + J£ as a function of Cp (top left panel), Cy (top right 
panel) and Ca (bottom panel). The cutoff scale has been taken to be Ay = lOOOGeF. 

(K(p')\s W b\B(p)) = , (18) 
where q = p — p' . Or alternately from the light cone sum rules [20] as: 

(K(p')\sj»b\B(p)) = {(p + pX - ^^ q ^m) + {^^q,}f P M 2 ) ■ (19) 
Note that we can relate these two sets of form factors by: 

f+ — f+ ) 

/- = (fo P ~ / + P ) • (20) 

In our numerical results we have followed the parameterization of Ball and Zwicky [20]: 

f = T2 
l-q 2 /m) it 
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J + l-g 2 /mf (l-g 2 /m 2 ) 2 
where the fitted parameters are given in Table 1. 





mi 




T2 


mfu 


f p + 


5.41 


0.1616 


0.1730 










0.3302 


5.41 



Table 1: The parameters for the B — > K form factors [19]. 



A. 2 The form factors for the B — > K* transition 

The form factors for the B — > K* transition can be written as [19]: 

(K*(p>)\s^b\B(p)) = 2V{q2) 



(K*(p')\s W b\B(p)) = t e ll (m B + m K *)A 1 (q 2 )-t(p + p') ll (e.q) 

-^(e.g)^[A 3 (g 2 )-A (g 2 )] , (22) 

where have again defined q = p — p'. For this transition we have used the parameterization of 
reference [19]: 

where 

^ = !M±I^ Ai(?2) _ , . (24) 

2mx* 2m^* 
Note that the fitted parameters used in the above equations have been given in Table 2. 
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